High-accuracy relativistic many-body calculations of van der Waals coefficients, Cq, 

for alkaline-earth atoms 
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Relativistic many-body calculations of van der Waals coefficients Ce for dimers correlating to 
two ground state alkaline-earth atoms at large internuclear separations are reported. The following 
values and uncertainties were determined : Ce = 214(3) for Be, 627(12) for Mg, 2221(15) for Ca, 
3170(196) for Sr, and 5160(74) for Ba in atomic units. 

PACS numbers: 31.10.-|-z, 34.20.Cf, 32.10.Dk, 31.15.Ar 



a. Introduction. The realization of Bose-Einstein 
condensation (BEG) in dilute ultracold samples of hy- 
drogen and alkali-metal atoms Li, Na, and Kb ||] 
has prompted a search for other atomic and molecu- 
lar species where BEG can possibly be attained. Of 
non-alkali atoms, so far the condensation was successful 
only with metastable helium ^. Gooling and trap- 
ping experiments with alkaline-earth atoms (in partic- 
ular Mg, Ga, and Sr) were recently reported (see, e.g., 
Ref. [§[ P, ||] ) and prospects of achieving the condensation 
with alkaline-earth atoms were also discussed 0. Ul- 
tracold alkaline-earth atoms possess several advantages 
over alkali-metal atoms. For example, utilization of the 
narrow spin-forbidden transition ^5*0 permits to 

optically cool atoms down to nano-Kelvin regime j|, ||. 
There is also a number of isotopes available with zero nu- 
clear spin, so that the resulting molecular potentials are 
not complicated by the "spaghetti" of hyperfine-structure 
states; this simplifies studies of trap losses and ultracold 
collisions 

Here we apply relativistic many-body methods to the 
determination of dispersion (van der Waals) coefficients 
Cq for the interaction of two identical alkaline-earth 
atoms in their respective ground states. The leading in- 
teraction of two ground-state atoms at large internuclear 
separations R is parameterized as —Cq/R^. Knowledge 
of the dispersion coefficients is required, for example, in 
determination of scattering lengths governing properties 



of Bose-Einstein condensates of dilute samples [10 . 

We employ several atomic relativistic many-body 
methods of varying accuracy. The dominant contribu- 
tion to Cg was evaluated with the configuration interac- 
tion (GI) method coupled with many-body perturbation 
theory (MBPT) |l^, |l^; smaller terms were computed 
using the less accurate relativistic random-phase approx- 
imation (RRPA) and Dirac-Hartree-Fock (DHF) meth- 
ods. The values were further adjusted with accurate the- 
oretical and experimental data for electric-dipole matrix 
elements and energies of principal transitions. We tabu- 



late the values of Ce for Be, Mg, Ga, Sr, and Ba. We also 
estimate uncertainties to be on the order of 1-2% for all 
alkaline-earth atoms, except for a 5% accuracy for Sr. 

b. Method of calculations. The dispersion coefficient 
Cq describes a second-order correction to energies of 
molecular terms related to the interaction of atomic 
dipoles induced by molecular fields. Due to large inter- 
nuclear separations involved, the two-center molecular- 
structure problem is reduced to the determination of 
atomic properties. The van der Waals coefficient may 
be expressed as Q 



Cq^6J2 



(E, - Eg) + [E, - Eg) 



(1) 



where and Eg are the wavefunction and energy of the 
atomic ground state, is an electric-dipole operator, 
and the summation is over intermediate atomic states 
'i'i and with respective energies Ei and Ej. Atomic 
units h = |e| = rUe = 1 are used throughout. The above 
relation can be recast into the Gasimir-Polder form 



[a{iuj)]'^ du! . 



(2) 



where a(«w) is the dynamic polarizability of imaginary 
argument defined as 



a{iuj) — 2Re 



{^g\D,\^,){^.,\D,\^g) 

{E, -^Eg)+iuj 



(3) 



The intermediate states in the sum, Eq. (|3|), can be 
separated into valence and core-excited states. We write 



a{iuj) = ay{iuj) + adiuj) 



j{iuj) . 



(4) 
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To determine the valence contribution ay we employ 
combined relativistic configuration interaction method 
and many-body perturbation theory (GI-I-MBPT). 
Smaller contributions of core-excited states ac are esti- 
mated with the relativistic random-phase approximation 
for the atomic core. In this method excitations of core 
electrons are allowed into the occupied valence shell and 
we introduce the correction acv to account for a sub- 
sequent violation of the Pauli exclusion principle; this 
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small correction is evaluated using the Dirac-Hartree- 
Fock method. 

Similar relativistic many-body techniques were in- 
volved in our previous high-precision determination of 
van der Waals coefficients for atoms with one valence 
electron outside a closed core Divalent atoms, 

considered here, present an additional challenge due to 
a strong Coulomb repulsion of the valence electrons. 
This strong interaction is treated here with the config- 
uration interaction method and smaller residual correc- 
tions (like core polarization) are treated with the many- 
body perturbation theory. The method, designated as 
CI+MBPT, was developed in Ref. [2[ |l|, ^ |l|, ^ . 

Here we briefly recap the main features of the 
CI-f-MBPT method. The complete functional space for 
electronic wavefunctions is partitioned in two parts: the 
model space spanning all possible excitations of the two 
valence electrons and an orthogonal space which adds 
various excitations of core electrons. The valence CI ba- 
sis set is saturated; e.g., the Ba ground state wavefunc- 
tion is represented as a combination of 1450 relativistic 
configurations in our calculations. Application of pertur- 
bation theory leads to effective operators encapsulating 
many-body effects and acting in the model space. For 
example, the CI valence wavefunctions are determined 
from the Schrodinger equation 



with the effective Hamiltonian defined as 
H,siE) = Ho + C + J:iE). 



(5) 



(6) 



Here Hq is the lowest-order Dirac-Fock Hamiltonian, C 
is the residual Coulomb interaction between valence elec- 
trons, and S is the energy-dependent self-energy opera- 
tor corresponding to core-polarization effects in model- 
potential approaches. The operator E completely ac- 
counts for the second order of perturbation theory. By 
the same virtue, one introduces an effective (dressed) 
electric-dipole operator Deff acting in the model space. 
We determine this effective operator using the random- 
phase approximation (RPA) ||2^, Qualitatively, the 
RPA describes a shielding of the externally applied field 
by the core electrons. 

The dynamic valence polarizability aylito) was com- 
puted with the Sternheimer |Q or Dalgarno-Lewis 
method implemented in the CI-I-MBPT-I-RPA frame- 
work. At the heart of the method is a solution of an 
inhomogeneous Schrodinger equation for a "perturbed" 
state \S^u;) 

iH,ff-Eg+iu;)\S-^^)^{D,U\^iJg), (7) 



so that 



a^iiu) = 2Re{^g\{D,%s\d-^^) . (8) 



In these expressions the electric-dipole operator Dcs is 
calculated at the CH-MBPT-I-RPA level of approxima- 
tion. Present approach is a frequency-dependent gener- 
alization of calculations of static dipole polarizabilities 



TABLE L Reduced matrix elements D and energy separations 
AEp for transitions from the lowest-energy nsnp^Pi to the 
ground ns^ ^5*0 state. 



D 



AEr, 





CI+MBPT+RPA 


accurate 


CI+MBPT 


Expt." 


Be 


3.26 


3.26(1)" 


0.194291 


0.193942 


Mg 


4.03 


4.03(2)" 


0.159173 


0.159705 


Ca 


4.93 


4.967(9)'' 


0.107776 


0.107768 


Sr 


5.31 


5.28(9)" 


0.098508 


0.098866 


Ba 


5.52 


5.466(23)" 


0.082891 


0.082289 



"Ref. |5|, 'This work, "Ref. '^Ref. fj. "Ref. 



reported in |2J]; technical details can be found in 
these works. 

The overwhelming contribution (on the order of 90%) 
to the value of the van der Waals coefficient, Eq. (|^), 
comes from the lowest-energy excited nsnp iPi° state. 
Therefore the calculated Cg are mostly sensitive to accu- 
racies of dipole matrix elements and energy separations 
of the principal nsnp ^Pf — ns^ ^Sq transitions. We ex- 
plicitly calculated these quantities using the same level 
of CI-j-MBPT+RPA approximation as employed in the 
solution of the inhomogeneous equation (^; these values 
are marked as CI+MBPT+RPA and CI+MBPT in Ta- 
ble |. We find a good agreement with more sophisticated 
ab initio ||2^, |l) and experimental values |Q, 25, |2^ 
(see TableHT) For Be we also computed additional many- 
body corrections; they can be neglected at the level of the 
quoted significant figures in Table |. We conservatively 
estimated an uncertainty in the matrix element for Be as 
a half of the difference between valence CI and correlated 
value. 

Due to the enhanced sensitivity of Cg to uncertainties 
in the dipole matrix element and the energy separation 
AEp of principal transitions, we further correct the calcu- 
lated dynamic polarizability by subtracting the ab initio 
CI+MBPT+RPA contribution of the principal transition 



ap{iu) = - 



AEr, 



3 (AEp) 



\{ns^^Sa\\D\\nsnp^P°)\^ (9) 



from a{iLu) and adding it back with experimental energies 
and high-accuracy matrix elements compiled in Table |. 

The "perturbed" state \S'i!i^) in Eq. (|^) is defined in the 
model space of the valence electrons, i.e., it is comprised 
from all possible valence excitations from the ground 
state l^'g). Since the core-excited states do not enter 
the model space, their contribution to the polarizability 
has to be added separately. Here we follow our work | ^ 
and use the relativistic random-phase approximation 1 2^ 
to determine the dynamic core polarizability as 



adiuj) 



to ('^m) 



(10) 



Here the summation is over particle-hole excitations from 
the ground state of the atomic core; are excitation en- 
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TABLE II: van der Waals coefficients C(, for dimers corre- 
lating to ground states of alkaline-earth atoms in a.u. Val- 
ues marked ab initio were determined in the relativistic 
CH-MBPT-I-RPA framework. The values marked final are 
ah initio values adjusted for accurate dipole matrix elements 
and energies of principal transitions, compiled in Table |. 



Be 



Me 



Ca 



Sr 



Ba 



Ah initio 213 631 2168 3240 5303 

Final 214(3) 627(12) 2221(15) 3170(196) 5160(74) 

Other works 

Stanton ^ 216 648 2042 3212 

S&C Hl 220 634 2785 

M&K |M1 208 618 2005 

A&Ch^l 254 2370 

Stwallef^] 683(35) 



ergies and are the corresponding electric-dipole oscilla- 
tor strengths. Accounting for core excitations is essential 
in our accurate calculations, especially for heavier atoms. 
For example, for Ba they contribute as much as 15% to 
the total value of Cg . 

The particle-hole excitations summed over in Eq. (|1£ 
include Pauli-principle violating excitations into the oc- 
cupied valence shell. We explicitly subtract their contri- 
bution; this small correction acv{i^) is computed with 
the Dirac-Hartree-Fock method. 

Our calculated dynamic polarizabilities satisfy two im- 
portant relations: (i) a{u) = 0) is the ground-state static 
dipole polarizability and (ii) as a consequence of the non- 
relativistic Thomas-Reiche-Kuhn sum rule, at large fre- 
quencies a{iuj) N, where N is the total number of 
atomic electrons. Indeed, for Ca we obtain a(0) = 160 
a.u., while the experimental value is 169(17) a.u. For 
Sr we obtain 199 a.u. which is in agreement with the mea- 
sured value of 186(15) a.u. And, finally, for Ba the 
computed static polarizability of 273 a.u. also compares 
well with experimental value |^ of 268(22) a.u. Simi- 
larly, at large w, in our calculations the product cj^ a(ia;) 
approaches 3.99 for Be, 11.9 for Mg, 19.71 for Ca, 37.1 for 
Sr, and 54.01 for Ba; these asymptotic relativistic values 
are slightly smaller than the exact nonrelativistic limits. 

c. Results and theoretical uncertainties. We com- 
bine various parts of the dynamic polarizability, Eq. (jj), 
and then obtain dispersion coefficients Cq with a quadra- 
ture, Eq. (H). The resulting values of van der Waals co- 
efficients are presented in Table ^ In this Table, val- 
ues marked ah initio were determined in the relativistic 
CI-I-MBPT-I-RPA framework. The values marked final 
are ab initio values adjusted for accurate dipole matrix 
elements and energies of principal transitions, compiled 
in Table |. 

Different classes of intermediate states in Eq. (||) con- 
tribute at drastically different levels to the total values of 
dispersion coefficients. For example, for Ca, the principal 
4s4p "^P^—As^ transition contributes 85% to the values 
of Cg, remaining valence- valence excitations contribute 
8%, core-excited states contribute 8% and the counter 



term acv modifies the final result only by -0.4%. To esti- 
mate dominant theoretical uncertainties we approximate 
Cg as 



3 /"^ g /"OO 

— / [ap{iuo)Yduj ^ / ap[iuj)ar{iuj)duj 

I" Jo Jo 



(11) 



Here ap is a contribution of the principal transition 
Eq. (^, and ar = a'^ + Oic is a contribution of the re- 
maining valence states (a'„ = ~ oip) and core-excited 
states. From a direct calculation for Ca wc find that this 
approximation recovers 99.3% of the Cg obtained from 
the full expression (^). Based on Eq. (|ll|) the sensitivity 
of Cg to uncertainties 5D in the matrix element D of the 
principal transition is 



5dC^ « (4C|^ + 2Cr) 



5D_ 
~D 



(12) 



To evaluate the sensitivity of Cg to uncertainties in the 
residual polarizability wc follow Rcf . p6| . In the second 
term of Eq. (11) a narrow function a^uj) is integrated 
with a relatively broad distribution aj.{iLo). Therefore we 
can approximate that 



ap{iuj)aj.{iLo)dhj ~ ar{0) / ap{iu!)duj (13) 
Jo 



and the sensitivity of Cg is 



pr 5ar{Q) 
' ar(0) 



(14) 



The uncertainty in the residual static polarizability 
(5ar(0) is a sum of uncertainties in the contributions of va- 
lence states beyond principal transition 5a'^ (0) and core- 
excited states (5ac(0). The RRPA static dipole core po- 
larizabilities for alkali-metal atoms are known |l^ to be 
in a 1% agreement with those deduced from semiempir- 
ical analysis of Rydberg spectra; we approximate that 
(5qc(0) ~ 0.01q!c(0). Further we estimate that 5a'^[Q) ~ 
6ap{0), i.e. the difference of the contributions of the prin- 
cipal transition to static polarizability calculated with 
CI-|-MBPT-f RPA and accurate values compiled in Ta- 
ble |. 

The error bars of the final values of dispersion coef- 
ficients in Table |l| were calculated by adding the un- 
certainties SnCe and Sa^Ce in quadrature. For all con- 
sidered alkaline-earth atoms the uncertainty in Cg in- 
duced by errors in matrix elements of principal transi- 
tion, 6dCq, dominates over (Sa^Cg. The estimated total 
uncertainties are in the order of 1-2% for all alkaline-earth 
atoms, except for Sr where the accuracy is 5%. Similar 
error analysis for alkali- metal atoms has proven to be 
reliable; for example, for Cs the predicted Cg — 6851(74) 
a.u. was found to be in agreement with a value |37]| of 
6890(35) a.u. deduced from an analysis of magnetic-field 
induced Feshbach resonances and photoassociation data. 
However, we emphasize that in the case of alkali-metals a 



4 



number of independent high-accuracy data was available 
for the dominant principal transitions ensuring reliability 
of derived dispersion coefficients. This is not the case for 
alkaline-earth atoms. In our present calculation we rely 
on the quoted uncertainties of accurate dipole matrix el- 
ements listed in Table |. 

A comparison with other theoretical and semiempiri- 
cal determinations is presented in Table There is a 
reasonable agreement among different approaches for Be 
and Mg; results for Ca are less consistent due to a more 
significant role of correlations and core-excited states. 
Coupled-cluster calculations by Stanton were most 
elaborate among theoretical treatments. We find a good 
agreement with his predictions. Unfortunately, most of 
the authors do not estimate uncertainties of their meth- 
ods. One of the exceptions is Ref. where sum rules 
and Pade-approximants were used to establish bounds on 
Cg. For Ca, they found 2740 < Cq < 2830 a.u. However, 
large uncertainties of underlying experimental data were 
not included in these bounds (see also Ref. |^); this 
explains a significant deviation of our prediction for Ca, 



Ce = 2221(15) a.u., from constraints of Ref. |3|. 

d. Conclusion. We carried out relativistic many- 
body calculations of van der Waals coefficients Cg for 
dimers correlating to two ground state alkaline-earth 
atoms at large internuclear separations. The values were 
adjusted with accurate theoretical and experimental data 
for the electric-dipole matrix elements and energies of 
the principal transitions. It is worth emphasizing that 
the dispersion coefficients depend sensitively on electric- 
dipole matrix elements of principal transitions. As more 
accurate data for the matrix elements become available, 
for example from photoassociation experiments with ul- 
tracold samples, the van der Waals coefficients can be 
constrained further within our many-body approach. 
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